Validation of techniques to mitigate copper surface contamination in CUORE 
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Abstract 

In this article we describe the background challenges for the CUORE experiment posed by surface contamination of inert detector 
materials such as copper, and present three techniques explored to mitigate these backgrounds. Using data from a dedicated test 
apparatus constructed to validate and compare these techniques we demonstrate that copper surface contamination levels better than 
10~ 7 - 10~ 8 Bq/cm 2 are achieved for 238 U and 232 Th. If these levels are reproduced in the final CUORE apparatus the projected 
90% C.L. upper limit on the number of background counts in the region of interest is 0.02-0.03 counts/keV/kg/y depending on the 
adopted mitigation technique. 
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1. Introduction 

Neutrinoless Double Beta Decay (J3/30v) is an area of vigor- 
ous experimental activity with potential for profound impact on 
modern questions in fundamental physics |UJ[2][3]. Observa- 
tion of this decay would immediately imply lepton number vio- 
lation and would establish the neutrino as a Majorana fermion. 
If neutrinos are indeed Majorana fermions, a measurement of 
the fifiOv rate would probe the absolute neutrino mass scale and 
possibly reveal the neutrino mass hierarchy. 

Sensitivity to the non-degenerate inverted hierarchy of neu- 
trino masses is a standard benchmark for next generation 
/3/30v searches which demands that very low background levels 
- of the order of few counts per ton per year in the region of in- 
terest - be achieved. All aspects of the experiment, for example 
selection of materials, machining and handling of components, 
and assembly procedures must be scrutinized for background 
control. Validation of effective control measures and quanti- 
fying the residual background is often nearly as challenging as 
the underlying experiment. This paper focuses on aspects of the 
background control and validation activities for the CUORE ex- 
periment |4] [5] |6] Q . Specifically, we present a study of three 
techniques explored to mitigate background from residual sur- 
face radioactivity on structural materials in the detector, partic- 
ularly copper. 

2. Overview of the CUORE detector 

The CUORE experiment, currently under construction un- 
derground at Laboratori Nazionali del Gran Sasso (LNGS), will 
search for f3/30v of 130 Te. The signature of this decay is a peak 
in the energy spectrum centered at the Q-value of the transition, 
at about 2528 keV [8,9, 10 1. The experimental goals include a 
background level of ~1(T 2 counts/keV/kg/y in an energy win- 
dow of ~100 keV around the Q-value, denoted the region of 
interest (ROI), and a high-precision measurement of the spec- 
trum in that region. 

The apparatus, shown in Figure [TJ will consist of a close- 
packed array of 988, 5x5x5 cm 3 cubic Te02 crystals, amount- 
ing to 206 kg of I30 Te. These will be cooled inside a cryo- 
stat to around 10 mK. At this temperature the crystals func- 
tion as highly sensitive calorimeters, converting the energy de- 
posited in their volume to a measurable temperature change. 
The bolometers will be arranged in a compact cylindrical ma- 
trix of 19 towers, each tower will contain 13 planes of four crys- 
tals. A copper skeleton will provide the mechanical structure to 
hold the crystals in each tower. The array will hang in vacuum 
inside a copper cylindrical vessel closed on the top and bottom 
with copper plates. The copper skeleton, denoted collectively 
as the copper holder, will not touch the crystals directly, instead 
PTFE standoffs will secure the crystals. Components made of 
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Figure 1: Set-up of the CUORE experiment: the 988 bolometers arranged in 
a 1 9 towers array, hanging in vacuum inside nested copper cylindrical vessels 
and provided with lead shields. On the right a detail of one CUORE tower. 

material other than copper or Te02 make up a small fraction of 
the detector. These components, denoted collectively as small 
parts, include the thermistors used to read out the bolometric 
signal, the silicon heaters, used to check for gain variations, the 
glue used to attach the thermistors to the crystals, the PTFE 
standoffs, and the readout wires. A complete description of the 
CUORE detector can be found in J6). 

3. Bolometer background from surface and bulk radioac- 
tivity 

In this Section we discuss aspects of bolometer behaviour 
that influence their susceptability to background from decays 
of surface radio-impurities and describe some discriminators 
to distinguish between surface and bulk contamination. In the 
context of this paper, bulk contamination represents unwanted 
impurities distributed uniformly throughout the volume of a 
material. This contamination arises from impurities present 
in the raw material, or introduced during manufacturing. On 
the other hand surface contamination refers to impurities from 
the environment that adhere to, or become embedded in the ex- 
posed surfaces of the material or to bulk impurities diffused into 
the surface. Since some impurities, for example long-lived de- 
cay fragments of 238 U or 232 Th, can become implanted inside 
the surface of the material the contamination is not strictly a 
surface distribution but rather may be distributed in a thin layer 
near the surface. 

Although bolometers make excellent calorimeters ifTTl fT2l . 
achieving FWHM energy resolution as high as 3.9 + 0.7 keV at 
2614.5 keV IfTTl . the phonon signal at the heart of the method 
shows very little dependence on either the identity of the parti- 
cle responsible for the underlying energy deposition or on the 
position of the event in the crystal. Electrons, gammas, a- 
particles and nuclear recoils depositing the same amount of en- 
ergy at any point in the bolometer produce virtually identical 
pulses lfT3lfT4l . The lack of particle identification and the ab- 
sence of a surface dead-layer on the crystals means that a's as 
well as y6's and y's form a background if they deposit an amount 
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Figure 2: Left plot: multiplicity = 1 spectra from Monte Carlo simulations of Pb contamination in the crystal bulk (solid-line histogram), spread on a very shallow 
depth (dotted-line histogram) and on a deeper depth (dashed-line histogram) of copper surface. Right plot: Monte Carlo simulated spectrum obtained summing up 
the energy released by two events in coincidence for a 2l0 Pb surface contamination on the crystals. 



of energy in the bolometer equivalent to the PfiOv transition en- 
ergy. 

Only a few decays in the chains of 238 U and 232 Th emit and 
y-particles with sufficient energy to mimic a fiflOv candidate, 
the most dangerous being the 2614.5 keV y-ray from 2()8 T1 — 
a daughter of 232 Th. On the other hand all a decays in the 238 U 
and 232 Th chains involve transition energies much higher than 
the I30 Te Q-value; if the entire transition energy is absorbed in 
one bolometer these do not contribute to the f3/30v background. 
However, this is not true if the a decay occurs on the surface 
of the detector or of the facing inert materials. In this case the 
transition energy is shared among more than one element and if 
the correct amount of energy is deposited in any crystal on the 
trajectory of the decay fragments, the pulse from that crystal 
could mimic a /3/30v candidate. Thus degraded ar's from surface 
impurities form a potentially serious background for bolome- 
ters. Studies performed on data acquired with bolometric ex- 
periments preceeding CUORE - like MiDBD ifTSl and Cuori- 
cino 1 16 1 - indicated surface contaminations as the most perni- 
cious source of background lfT71[T8l . 

A very useful parameter to identify and veto decays from 
surface contamination is event multiplicity. This is defined as 
the number of crystals firing in coincidence, typically within 
100 ms of each other (typical rise and decay times are of the or- 
der of tents and hundreds of ms respectively, and the expected 
average rate is around a few mHz). Particles emitted from de- 
cay of surface contamination can hit more than one bolome- 
ter, whereas a true /3/30v candidate is expected to be contained 
within a single crystal with an efficiency of ~83% lfl9l . As an 
example, we show in Figure[2]and discuss below, results from a 
Monte Carlo (MC) simulation which illustrate the multiplicity 
characteristics of a decay of 210 Po, a daughter of 210 Pb. 

First we consider if the 210 Pb is present in the bulk or in the 
surface of Te02 crystals. In the bulk, the energy of both the a 
and recoiling nucleus are absorbed fully in the crystal. Thus a 
gaussian peak centered at the isotope transition energy in the 
multiplicity =1 (Ml) spectrum is found, see the solid-line his- 
togram in the left plot of Figure [2] 

On the other hand, in the surface region of a crystal, ar's - 
or/and the recoiling nuclei - are not fully stopped in one crystal 
and sometimes enter a second bolometer causing both to fire 



in coincidence. The case where two bolometers fire in coinci- 
dence corresponds to multiplicity=2 (M2) events. The summed 
energy of the two events is denoted M2Sum. The M2Sum spec- 
trum has a gaussian peak centered at the full isotope transition 
energy since the a and recoiling nucleus are fully absorbed in 
the two crystals. See right plot of Figure[2] For more details see 
Ref. 11201 . 

Next we consider the case of contamination in the surface 
region of inert materials such as copper. Degraded or's may es- 
cape the copper and hit a bolometer. Only the energy deposited 
in the bolometer can be observed, the energy deposited in the 
inert material is lost. For very shallow surface contamination 
a peak at the a energy — the transition energy minus the nu- 
clear recoil energy — is visible in Ml spectrum, see the dotted- 
line histogram in the left plot of Figure [2] However, when the 
contamination is deeper, the range of energies possible for the 
degraded a becomes broader and an almost flat continuum, ex- 
tending from the full a energy down to the lowest energies, is 
found. For example, see the dashed-line histogram in the left 
plot of Figure [2] 

The energy interval from 2.7-3.9 MeV is an ideal region 
to study the degraded a background as it is above the high- 
est y-line in the 238 U and 232 Th chains, 2614.5 keV, but be- 
low the lowest cr-line at about 3947 MeV from 232 Th. The 
only peak visible in Te02 experimental spectra in this region 
is at 3270 keV, this line comes from 190 Pt bulk contamination 
of Te02 crystals due to inclusion of platinum atoms from the 
crucible during crystal growth BTTl . In the subsequent Sections 
we refer to the energy interval from 2.7 MeV to 3.9 MeV, ex- 
cluding the 200 keV window centered on the 190 Pt peak, as the 
degraded a window. 

The Cuoricino experiment |[T6l[T8l . a precursor to CUORE, 
found that the observed background in the ROI for the /3/30v of 
130 Te, 0.18 + 0.01 counts/keV/kg/y, was consistent with the 
expectation from three classes of sources: (i) multi-Compton 
events caused by 2614.5 keV y-radiation from 208 T1 supported 
by 232 Th contamination of the cryostat or its shields; (ii) ra- 
diation from surface contamination of the Te02 crystals with 
232 Th, 238 U, or 210 PbQ and (iii) radiation from surface contam- 
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ination of inert materials surrounding the crystals, most likely 
copper. In the case (ii) and (iii) the radiation is primarily de- 
graded a's with a small contribution from (3+y events. Using 
signatures from outside the ROI to fix the normalization of each 
source the extrapolated count rate in the ROI was found equiv- 
alent to: (i) 30 + 10%, (ii) 10 + 5% and (iii) 50 + 20% of the 
observed background rate. Although this finding does not ex- 
clude other sub-dominant backgrounds, the study indicates that 
~30% of the background originated from the cryogenic appara- 
tus and ~70% from detector surface contamination, particularly 
degraded a's. 

To meet the ROI-background goal of CUORE both these 
sources must be addressed. Materials for the CUORE cryo- 
stat and the internal and external shields with sufficient radiop- 
urity to ensure a ROI-background contribution far below the 
target level of 10~ 2 counts/keV/kg/y have already been iden- 
tified Jfj). For the Te02 crystals and the copper holder the ac- 
ceptable 238 U and 232 Th contamination levels are as low as 1 
/iBq/kg for the bulk and 1-10 nBq/cm 2 for the surface. While 
copper bulk contamination can be adequately validated with 
HPGe spectroscopy and Neutron Activation Analysis (NAA), 
these techniques do not have sufficient sensitivity to validate 
the radiopurity of the copper surface, the Te02 surface, and the 
Te0 2 bulk. 

To meet this validation challenge a series of test arrays con- 
taining a few bolometers arranged in the style of one or a few 
floors of CUORE were operated in dedicated runs at LNGS Q. 
The results of Te02 crystal surface validation runs were re- 
ported in (20); 90% C.L. upper limits on 238 U and 232 Th surface 
concentration corresponding to 3.8 nBq/cm 2 and 2.0 nBq/cm 2 
respectively were demonstrated. The role of copper and the re- 
sults and interpretation of test runs dedicated to copper surface 
validation are discussed in the remainder of this paper. We fo- 
cus exclusively on copper since next to Te02 this makes up the 
largest material surface in the detector. The surface area of the 
detector small parts is so small that in most cases their contribu- 
tion to the ROI-background is already constrained by Si barrier 
detector measurements to be a negligible fraction of the copper 
contribution. 

4. Copper contamination and treatment for CUORE 

The mechanical and cryogenic properties of copper make it 
an ideal material for CUORE. NOSV copper [22], produced by 
the Norddeutsche Affinerie |f23| , was selected for the experi- 
ment because of its high residual-resistivity ratio (RRR) (certi- 
fied to be higher than 400) — a constraint from our cryogenic 
application — and because samples were found to be extremely 
radio-pure. HPGe and NAA measurements limit bulk contami- 
nation to less than 5 pg/g of 238 U and 0.5 pg/g of 232 Th l24ll . 

Cosmogenic activation of the copper, especially 60 Co pro- 
duction, can be controlled to an acceptable level by storing it 
underground at LNGS except for the time needed for machin- 
ing and cleaning. We estimate the time above ground of about 4 



naturally occuring isotopes: Th, U. Pb is also studied as an indepen- 
dent source as it is often found out of secular equilibrium. 




(a) (b) 

Figure 3 : Photograph of the Three Towers Test detector with copper box (a) and 
picture of the detector without copper coverage (b). From top to bottom: the 
polyethylene tower (Tl), the LNGS tower (T2), and the Legnaro tower (T3). 



months from cast to final storage. With this exposure time, us- 
ing the production cross sections in [25 1, we expect the cosmo- 
genic activity of the copper holder will be less than 50 juBq/kg 
when CUORE starts taking data, ensuring a negligible contri- 
bution to the background in the /5j60v-ROI. 

Although the copper bulk-contamination satisfies the re- 
quirements of CUORE the surface contamination of commer- 
cially available copper needs to be reduced. Unfortunately stan- 
dard surface cleaning procedures, generally based on pickling 
and etching, appear to be unable to reach the desired radiopu- 
rity. For Cuoricino the surfaces of all the copper components 
facing the bolometers were treated at Laboratori Nazionali di 
Legnaro (LNL) with a procedure similar to one used for res- 
onant cavity production |26l [27). It involved a sequence of 
tumbling, to reduce the surface roughness; chemical etching, 
to reduce chemical contamination; and passivation, to decrease 
the possibility of re-contamination. Although this procedure 
proved more effective than standard methods lfT31[T7l . the rate 
in the ROI due to contamination of copper surfaces in Cuori- 
cino was constrained to be 50 ± 20% of the observed ROI- 
background rate|^] 

To improve on this for CUORE we strive to both mini- 
mize the amount of copper and other inert materials facing the 
bolometers and to identify better surface treatment to mitigate 
the background from the remaining surfaces. Three surface 
treatments techniques were chosen to test: (i) wrapping of sur- 
faces with polyethylene, (ii) simple surface cleaning with ultra- 
clean acids, and (iii) a modification of the LNL procedure. A 
dedicated bolometric test, called the Three Towers Test, was 
organized to validate and compare these treatments. 



"Important environmental muon and neutron contribuions were ruled-out on 
the basis of simulations and measurements [28 17 29 1 
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Figure 4: Photographs of the three 12-crystal arrays (the three towers ) without 
their cylindrical shields: (a) the polyethylene tower (Tl), (b) the LNGS tower 
(T2), and (c) the Legnaro tower (T3). 

5. The Three Tower Test (TTT) 

The apparatus, known as the TTT detector, consisted of three 
12-crystal arrays separated from each other by copper shields, 
see Figures[3]and[4] The 36 crystals were all taken from Cuori- 
cino production. Before being installed in the TTT detector 
their surfaces were retreated using new procedures developed 
for CUORE EUll . Great care was taken to ensure the crystals 
had similar contamination levels, history, and treatment. The 
copper of the three towers was taken from one single batch 
of NOSV copper produced by the Norddeutsche Affinerie, and 
machined following identical procedures. The surfaces of the 
copper parts were treated with one of the three following pro- 
cedures □ 

• Tl: For the first or top tower, all the copper pieces were 
simply cleaned with soap, treated with a mixture of H2O2, 
H2O and citric acid, and then wrapped with a few layers, 
~70 |/m thick on average, of polyethylene film. Given that 
the range of a particles in polyethylene is about 20 yum, 
this is sufficiently thick to stop any a emission from the 
copper surface. Commercial polyethylene film was cho- 
sen because it is relatively easy to handle; measurements 
of samples showed good radiopurity, <100 ppt of 238 U and 
<700 ppt of 232 Th; and the procedure achieved good re- 
sults in previous tests. 

• T2: For the second or middle tower, a purely chemical pro- 
cess with ultra-pure reagents was used. The first step was 
simple cleaning with soap and water. The copper pieces 
were then subjected to electroerosion with 85% phospho- 
ric acid, 5% butanol, and 10% water, followed by chemical 
etching with hyper-pure nitric acid. Finally a passivation 
step with H2O2, H2O and citric acid was performed. The 
procedure was developed at LNGS on the basis of 1 30 1 . 

• T3: For the third or bottom tower, the LNL cleaning pro- 
cedure was modified to include: Tumbling, Electropol- 



3 The specifications of the reagents used in the threatments are: soap - Micro 
90, Sigma Aldrich [31]; phosphoric acid - 85%, RPE analytical grade, Carlo 
Erba (32); H 2 - obtained with Milli-Q system (33], 18.2 Mfl*cm at 25°C; 
solid citric acid - 99%, Sigma Aldrich 1311 : nitric acid - 67-69% RS hyper- 
pure distilled with DuoPur purification system from Milestone, Carlo Erba 1 32]; 
Butanol - RPE analytical grade, Carlo Erba 1 32]. 
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3.0 ± 0.5 


23 ±5 


T2 


1 




4.5 ± 1.2 


12 ±3 


T3 


1 




5.2 ± 1.0 


11 ±3 


Tl 


2 


2.6+ 1.0 


6.2+ 1.0 


21+4 


T2 


2 


1.6 ±0.3 


3.7 ± 0.7 


27 ±5 


T3 


2 


3.6+1.0 


4.2 ± 1.5 


17 ±4 



Table 1 : Average FWHM measured by the Three Towers in the two data sets. 

ishing, Chemical etching and Magnetron plasma etch- 
ing (TECM). Final extraction and the following storage 
of parts was performed in a clean room to avoid re- 
contamination of the magnetron treated surfaces. Ultra 
clean reagents were not employed since the expense of 
supplying such a large TECM plant was deemed unsus- 
tainable and the procedure is designed for zero-deposition 
of foreign material on the treated piece. 

The TTT apparatus was installed in the cryostat used for 
Cuoricino [18|, the setup was identical to Cuoricino with re- 
spect to cryogenics, shields and front-end electronics, see ref- 
erences lfT~8l and Ifl9l . The data acquisition system (DAQ) de- 
signed for CUORE was used. The DAQ digitizes pulses from 
each bolometer front-end channel with an 18-bit ADC sampling 
at 125 Hz. The system acquires 626 samples, corresponding to 
5.008 s, for each triggered pulse. Data were collected in two 
campaigns: data-set 1 (Dsl) which ran from September to Oc- 
tober 2009 and data-set 2 (Ds2) which ran from October 2009 
to mid January 2010. During Dsl a 40 K source was present to 
validate parts of the offline analysis, this source was removed 
for Ds2. The detector was operated at a higher energy thresh- 
old during Dsl to reduce the trigger rate from the 40 K source to 
an acceptable level. 

Unfortunately not all the electrical connections survived the 
apparatus cool down, ultimately 12, 7, and 7 crystals were read- 
able in Tl, T2 and T3 respectively. The readable detectors 
performed quite well, the average FWFLM resolution measured 
from the total spectrum (i.e. no multiplicity cuts applied) in 
each tower is shown in Table Q] for both data sets. Three en- 
ergy regions are considered: (i) low energy, i.e. the 352 keV 
214 Pb peak; (ii) an intermediate region near the ROI, i.e. the 
2614.5 keV 208 T1 peak; and (iii) the a region, i.e. the 5407 keV 
210 Po peak. The first region was below threshold for Dsl. The 
poorer average resolution which is evident at higher energy, the 
a region, is not due to the resolution of the individual bolome- 
ters but rather to the inter-bolometer calibration within a sin- 
gle tower. The highest available calibration peak is in fact 
2614.5 keV, beyond this point the calibration must be extrap- 
olated. 

5.7. Data analysis 

Data from the TTT were analysed using the software frame- 
work developed for CUORE and already successfully used for 
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Cuoricino |19|. The off-line analysis identifies and rejects oc- 
casional periods with excess noise, applies pulse-shape analy- 
sis to identify and remove spurious pulses, and evaluates pulse 
height. Pulse height is a proxy for energy, the precise relation- 
ship between pulse height and energy is fixed using calibration 
sources. The off-line analysis also identifies time coincidences 
between events for multiplicity analysis. Ultimately the multi- 
plicity analysis was only effective for Tl during Ds2 because 
this was the only tower where all 12 crystals were readable and 
the event rate from the 40 K source was too high during Dsl. 

The algorithms used for TTT analysis were identical to those 
described in Ifl9l except a new method used for detector re- 
sponse stabilization Il34ll . Detector response stabilization is a 
correction applied to pulse amplitudes to account for gain in- 
stabilities and temperature fluctuations. The new algorithm was 
validated using the 1461 keV y-peak from the 4l) K source data 
in Dsl. 

For each tower the average efficiencies are evaluated inde- 
pendently for the two data-sets following the procedure de- 
scribed in [ 19 1 . All the rates discussed in the following Sections 
are corrected for the efficiency. 

Table [2] summarizes the number of working detectors, the 
massxtime exposure, the thresholds and the possibility of ap- 
plying a multiplicity analysis for each data set. 

Next we introduce four classes of energy spectra which are 
important for the data analysis. 

• Total spectrum: Each event is recorded at the correspond- 
ing energy in the spectrum regardless of any coincidence 
between crystals. 

• Ml spectrum: Events are selected if no other event oc- 
curred in the other bolometers in the tower within the co- 
incidence window. A multiplicity =1 is assigned to such 
events. This is also referred to as the anti-coincidence 
spectrum. 

• M2- and M2Sum-spectrum An event is selected if there 
is exactly one more event in one of the other bolometers 
in the tower within the coincidence window. A multiplic- 
ity=2 is assigned to the two events and they are said to 
satisfy the coincidence cut. In the M2 spectrum the in- 
dividual energies are recorded separately. In the M2Sum 
spectrum an entry is recorded at the energy corresponding 
to their sum. 

• M>2: Multiple coincidence cut. An event is selected if 
there are two or more events in the other bolometers in the 
tower within the coincidence window. 

Interpretation of the data in terms of higher level quantities 
such as the surface distribution of impurities relies heavily on 
Monte Carlo simulation. We use the Geant4 toolkit [ 35 1 with 
our detector geometry implemented to track particles, their in- 
teractions and energy deposition. Event generation relies on 
a package called GENDEC, developed and tested for previous 
TeC»2 bolometric arrays [ 17], to generate particles with the cor- 
rect time correlation along the decay chains of isotopes of in- 
terest. Another package called G2TAS [17| simulates specific 



Working detectors 



Tower 


Dsl 


Ds2 


Tl 


12 


12 


T2 


7 


7 


T3 


7 


7 


Exposure [kgxy] 


Tower 


Dsl 


Ds2 


Tl 


0.87 


1.09 


T2 


0.51 


0.76 


T3 


0.52 


0.74 


Threshold 


keV] 


Tower 


Dsl 


Ds2 


all 


500 


100 


Multiplicity 


Tower 


Dsl 


Ds2 


all 


NO 


only Tl 



Table 2: Number of working detectors, exposure, energy threshold and enabling 
of multiplicity analysis for the three towers in the two data sets (Dsl and Ds2). 



Tower 


Th 


U 




238 keV peak 


352 keV peak 




[c/h]xl0- 3 


[c/h]xl0~ 3 


Tl 


5.6 + 2.3 


5.4 ±2.1 


T2 


9.5 ± 3.5 


6.6 ± 2.7 


T3 


5.4 ± 2.5 


3.1 ± 1.6 



Table 3: Intensity of the 238 keV and 352 keV peaks (due to 232 Th and 238 U 
chains respectively) in the three towers. No coincidence cuts have been applied. 
Only Ds2 has been used. 

acquisition and analysis features of bolometers arrays such as 
pile-up, event multiplicity, and detector thresholds and resolu- 
tions. The complete simulation chain used in this work is identi- 
cal to that described in detail in ifTTIl . except that we adopted the 
recently developed Livermore Physics Lists 13 61 . which better 
implements low energy interactions. In IfTTl it was shown that 
calibration measurements obtained by exposing a Te02 bolo- 
metric array to a 232 Th wire source were reproduced at better 
than 5% with MC simulations. Given that the uncertainties in 
the Geant4 physics models are estimated to be about 5% [37| 
and particle transport, interaction, and detection add another 
few percent, on the whole we consider 10% as a conservative 
estimate for systematic uncertainty in the simulation output. 

5.7.7. Comparison of surface radiation signatures measured in 
each tower 

Three signatures are studied to constrain surface contamina- 
tion with 238 U, 232 Th, and 210 Pb. Namely the degraded a sig- 
nature between 2.7 MeV and 3.9 MeV described in Section [3] 
the 238 keV y-line from 232 Th chain and the 352 keV y-line 
from 238 U chain. These low energy gammas have a short range 
(about 1 cm for a 300 keV gamma in copper) and thus events 
from these lines can only come from bulk and surface contam- 
ination of the crystals or from bulk and surface contamination 
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Figure 5: Comparison of the TTT spectra in the a region. Tl spectrum is shown in the top plot, T2 in the middle and T3 in the bottom one. No coincidence cuts are 
applied. The complete population of events has been used. 



Tower 


continuum 




2.7-3.9 MeV 


Tl 


0.068 ± 0.006 


T2 


0.120 ±0.012 


T3 


0.072 ± 0.008 



Table 4: Counting rates [counts/keV/kg/y] recorded in the 2.7-3.9 MeV energy region with the three towers. Data refer to the sum of the two data sets, no coincidence 
cuts are applied. 



of inert materials in their proximity. They do not have sufficient 
energy to reach the detectors from far inert materials. 

Table [3] shows the measured rates for the low energy y-lines 
for Ds2; both lines were below threshold in Dsl. Based on 
HPGe measurements of copper and bolometric measurements 
of the crystals we estimate that less than 25% comes from bulk 
contamination with 238 U and 232 Th (see Section [4] and ||20l ). 
Therefore the majority of these events are due to surface radi- 
ation from the crystals or copper. The measured intensities in 
each of the towers are consistent within the uncertainty. 

Figure [5] shows the Total spectrum summed over Dsl and 
Ds2 for Tl, T2 and T3. The rates in the degraded a window 
are reported in Table [4] The Tl- and T3-rates are compatible 
within lcr, whereas the rate in T2 is 70% higher. To benchmark 
and compare the effectiveness of the copper surface treatments 
the relative contributions of crystal surface contamination vs. 
copper surface contamination to the observed Total rates must 
be quantified. To do this directly requires a multiplicity anal- 
ysis; this is done for Ds2 of Tl in the next Section where we 
demonstrate that less than 20% of the observed continuum rate 
is due to surface radiation from the crystals. Multiplicity analy- 
sis was not possible for T2 and T3 due to loss of channels dur- 
ing cooldown. Instead we rely on the fact that all 36 crystals, 
all belonging to the same production batch and treated identi- 
cally, should have similar bulk and surface contamination lev- 
els. Thus while crystal surface contamination may contribute 
to the continuum counting rate, its variation between the tow- 
ers should be dominated by differences in the residual surface 
radiation from the treated copper surfaces. 



5.1.2. Constraints on the contribution of crystal surface radia- 
tion to the 2.7-3.9 MeV rate for Tl 

The Ml, M2, M2Sum and M>2 spectra used for the multi- 
plicity analysis are shown in Figure [6] 

In the Ml spectrum we note the following features: the 190 Pt 
peak at 3270 keV, peaks between 4 and 5 MeV identified with a 
transistions in the chains of 238 U and 232 Th, and peaks at about 
5407 keV from a decay of 2I0 Po (a daughter of 2I0 Pb). The 
only peak centered at the energy of the a fragment of any of the 
a transitions considered is the peak at 5304 keV. This is con- 
sistent with a very thin layer (of the order of 0.01 fim) of 21() Pb 
contamination on the surface of the crystals or inert materials 
facing them such that the a is detected in a facing bolometer 
without any loss of energy while the recoiling nucleus remains 
in the source element. The clear rate discrepancy between the 
intensity of the 210 Po a peaks and those identified with the 238 U 
chain indicates that the 210 Pb contamination is not primarily 
supported by 238 U. A break in the secular equilibrium in the 
lower part of both the 238 U and 232 Th chains is also evident 
from the discrepancy between the intensity of the peaks above 
6 MeV and those in the 4-5 MeV range (see also rates in Table 

In the M2 spectrum the only visible structure is at 5304 keV. 
It is an asymmetric peak, with a tail on the left side. This is 
consistent with 210 Po decays occurring in a shallow layer on 
the crystal surface where the a particle and recoiling nucleus 
release energy in two facing detectors. 

In the M2Sum spectrum the most clear peak is at 5407 keV. 
This peak is a different signature of the same 21() Pb contami- 
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Figure 6: Ml (at the top), M2 (the second from the top), M>2 (the third from the top) and M2Sum (the bottom one) energy spectra in the a region for Tl acquired 
events. The name of the parent of the decay chain is written above highest peaks. 



Mult. 


continuum 


U/Th 


Po/Pb 


U/Th 




2.7-3.9 MeV 


4-5 MeV 


5-6 MeV 


6-8 MeV 


Total 


0.068 ± 0.009 


0.29 ± 0.02 


1.38 ±0.07 


0.03 ± 0.005 


Ml 


0.052 ± 0.008 


0.28 ± 0.02 


1.30 ±0.07 


0.02 ± 0.004 


M2 


0.009 ± 0.003 


0.0018 ± 0.005 


0.09 ± 0.01 


0.005 ± 0.002 


M>2 


0.008 ± 0.003 


0.001 ± 0.001 


0.008 ± 0.003 


0.005 ± 0.002 



Table 5: Rates [counts/keV/kg/y] measured in Tl in the a region with different multiplicity constraints. Only the statistics collected in the second part of the 
measurement is used. 



nation on the crystal surface responsible for the 5304 keV line 
in the M2 spectrum. Using a MC simulation we extrapolate 
the intensity of the M2Sum peak at 5407 keV to estimate the 
contribution of crystal surface contamination to the Ml peak at 
5304 keV. The extrapolated rate is about a factor 5 lower than 
the measured one. Thus most of the events in the 5304 keV 
peak of the Ml spectrum are from 210 Pb contamination at a 
shallow depth in the inert material facing the crystals, mainly 
polyethylene in the case of Tl. 

In the M>2 spectrum no peaks are evident. Only a few 
sources are expected to produce high-multiplicity multiplets of 
events with a continuum energy spectrum. For example muon 
cascades, neutron interactions, y, and [5+y from decays with 
high Q-values (like 208 T1 decay whose Q-value is of about 
5 MeV) taking place in the detector proximity. 

Table [5] lists the rates measured in the degraded-ff region in 
Ds2 after the different multiplicity cuts. In addition, the re- 
sults for three higher energy intervals are shown. These in- 
tervals cover transition energies of a decays in the chains of 
238 U, 232 Th (4-5 MeV and 6-8 MeV) and 210 Pb (5-6 MeV) and 
their daughters. Thus events in these intervals give a first indi- 
cation of the presence of a contamination in the crystals bulk 
and/or surface. As explained in Section [3] a contribution to the 
peaks at the a energy in the Ml spectrum could arise also from 



surface contamination distributed in a very shallow depth of the 
facing materials, but this is expected to give a minor contribu- 
tion. 

Following the methodology in |20|, we use simulation to es- 
timate the maximum contribution to the 2.7-3.9 MeV region 
from 238 U, 232 Th and 210 Pb impurities on the crystal surfaces. 
For each chain, the simulation is used to predict the Total, Ml, 
M2 and M2Sum spectra for an assumed surface distribution of 
the impurity. The surface impurity distribution is taken of the 
following form: p = poe~ x ^, where po is the contamination 
per unit volume at the surface, x is the distance into the crys- 
tal, and A is the mean penetration depth of the impurity. For a 
given A the comparison between MC simulation and measure- 
ment rates fixes the normalizations of the Total, Ml, M2 and 
M2Sum simulated spectra for each source. The integral of p 
over the volume, multiplied for the evaluated normalization fac- 
tor and divided by the surface of the source element (on which 
the distribution is assumed to be uniform) gives what we quote 
as surface contamination 8, i.e. the concentration of impurities 
measured in Bq/cm 2 . 

For the 238 U and 232 Th chain, with A fixed, we estimate the 
maximum 6 consistent with the rates measured in the Ml- and 
M2Sum-spectra at the Q-values of the a transitions in the chain. 
Where no peaks were visible in the spectrum, the 90% C.L. 
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Feldman-Cousins upper limit 1 38 1 is evaluated using the num- 
ber of counts recorded in a 60 keV energy interval centered 
at the transition Q-value and assuming a background of zero 
counts We repeat this procedure for a range of values of 
A. A similar procedure is used for 210 Pb, except that the es- 
timate for 6 is only required to be consistent with the rate 
measured at the 5.4 MeV peak in the M2Sum spectrum. De- 
pending on the choice of A, we find values for 6 ranging from 
3xl(T 9 Bq/cm 2 to 9xl(T 9 Bq/cm 2 for 238 U and 232 Th, and from 
4 xl(T 8 Bq/cm 2 to 8xl(T 8 Bq/cm 2 for 210 Pb. 

Given a pair (A, 6), the expected count rate in the degraded a 
window in the Ml or Total spectrum is fixed by the simulation. 
The maximum expected count rate from 238 U, 232 Th and 210 Pb 
combined ranges between 2% and 10% of the measured rate 
for the Ml spectrum and between 3% and 20% of the measured 
rate for the Total spectrum. Again the range corresponds to the 
interval explored for A, with the lowest percentage correspond- 
ing to the smallest A. We conclude that less than 20% of the 
rate measured between 2.7-3.9 MeV in the Total spectra for the 
towers (see Table |4| is from crystal surface contamination. 

For two reasons we neglect the possible contribution of en- 
vironmental neutrons and muons to the difference in mea- 
sured rates. Firstly neutron and muon backgrounds should be 
the same for all three towers since they are in the same lo- 
cal environment. Secondly, on the basis of data from previ- 
ous TeC>2 experiments the expected contributions to the de- 
graded a window of the Total spectrum are negligibly low: 
the 90% C.L. upper limit for the expected neutron contribution 
is 8xl0~ 6 counts/keV/kg/y ifTTl and the expected contribution 
from muons is (1.26 ± 0.3)xlO~ 3 counts/keV/kg/y J29J. The re- 
maining candidate source is surface contamination of inert ma- 
terials surrounding the crystals, mainly the copper holder for T2 
and T3 and the multiple layers of polyethylene film for Tl. 

Since the rates measured in Tl and T3 are comparable we 
conclude that the polyethelyene and TECM treatments are sim- 
ilarly effective. The treatment for T2, which had a 70% higher 
rate than the other two is the least effective. 

5.7.5. Copper surface contamination derived from the mea- 
sured Total rates 

In this Section we proceed with the evaluation of the surface 
contamination of the copper facing the detectors. We do this 
only for T2 (chemical etching) and T3 (TECM) where the cop- 
per faces the crystals without any interposed material. 

We focus on the following signatures to estimate the concen- 
tration of 238 U, 232 Th and 210 Pb impurities in the surface of the 
copper holder from the measured spectra: (i) counting rate in 
the degraded a window, (ii) the 238 keV and 352 keV peak in- 
tensities, (iii) the intensities of a peaks centered at the a energy 
(transition energy - nuclear recoil energy). None of these signa- 
tures is exclusively due to contamination of the copper surface 
(for example impurities of the crystals themselves could con- 
tribute) and in particular, in the case of the degraded a window, 



The choice of considering background counts leads to an overestimation 
of the crystal contamination and therefore to a more conservative evaluation. 



Source 


T2 [Bq/cm 2 ] 


T3 [Bq/cm 2 ] 


232 Th 


1.3x10-' 


6.8xl0~ 8 


238 u 


1.3x10-' 


6.5xl0- 8 


21Up b 


8.8xl0~ 7 


8.6xl0~ 7 



Table 6: Evaluated 90% C.L. upper limits for T2 and T3 copper surface con- 
tamination. 

more than one species can contribute. For a quantitative eval- 
uation we assume that the 100% of the measured rate is due to 
each species of impurity on the copper surface, 238 U, 232 Th and 
210 Pb, in turn. The surface concentrations thus derived are con- 
servative upper limits. The main steps of the analysis are the 
following: 

• For each species we simulate events assuming an expo- 
nential depth profile for the contamination as described in 
the previous Section. We consider different characteristic 
penetration depths, A, ranging from 0.1 pm to 10 pm. The 
former is the most shallow distribution whose effects can 
be experimentally identified and the latter is the range of a 
5 MeV a in copper. Impurities at greater depths are indis- 
tinguishable from bulk contamination. 

• For each of the considered signatures we evaluate the de- 
tection efficiency, i.e. the fraction of simulated decays that 
result in an event with the signature of interest. We then 
identify the most efficient signature of the contaminant be- 
ing studied. 

• For the most efficient signature the maximum activity 
compatible at 90%C.L with the experimentally measured 
counting rate is calculated. The copper activity per unit 
area is then calculated by dividing this result by the emit- 
ting surface area. 

Since the actual characteristic depth is unknown, we report in 
Table [6] the contamination upper limits corresponding to the 
depth profile which gave the worst, i.e. the highest, value. The 
meaning of this upper limit is that in the worst configuration the 
copper surface contamination is below the reported activity at 
the 90% C.L. 

6. Contribution of copper to the ROI for CUORE 

An important final step is to estimate the residual background 
rate in the ROI for CUORE to confirm if it is compatible with 
the goal of the experiment (i.e. 10" 2 counts/keV/kg/y). As- 
suming a completely functioning array, which will allow full 
exploitation of multiplicity analysis, the figure of interest for 
CUORE is the background in the ROI of the anti-coincidence 
or Ml -spectrum. 

For a given tower of TTT the main observable used to bench- 
mark the treatment is the rate in the degraded a window (2.7 
to 3.9 MeV) of the Total spectrum. To estimate the transfer 
coefficient that connects this observable to the the Ml -rate of 
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the CUORE spectrum in the ROI we turn again to MC simula- 
tion. For a given distribution of impurities on the surface of the 
copper structure we simulate datasets for both the TTT tower 
and CUORE and extract both rates of interest, i.e. the rate in 
the 2.7-3.9 MeV window of the Total spectrum for TTT and the 
rate in the /?/30v-ROI of the Ml spectrum for CUORE. The ratio 
of these rates is an estimator for the transfer coefficient for the 
associated surface distribution. Due to the completness of the 
simulation this estimator simultaneously includes: the geomet- 
rical effects, i.e the fact that on average less copper faces the 
crystals in CUORE than in TTT; the anti-coincidence effects, 
i.e the fact that CUORE will benefit from full exploitation of 
multiplicity analysis; and spectral effects, i.e the fact that the 
degraded alpha spectrum increases approximately linearly with 
energy and thus has a lower contribution in the f3/30v-ROl rel- 
ative to the 2.7-3.9 MeV window. Given the transfer coeffi- 
cient the associated ROTbackground in CUORE can be esti- 
mated from the experimentally measured rate in the degraded 
alpha window of the TTT tower. Transfer coefficients were 
estimated separately for 238 U, 232 Thand 210 Pb surface distribu- 
tions and in each case, as before, a range of values for the mean 
penetration depth A are used. From this set of transfer coeffi- 
cients we choose the one that gives the worst ROI-background 
for CUORE . 

Following this procedure we estimate the ROI background 
in CUORE for the polyethelyene treatment, Tl, and the TECM 
treatment, T3. For Tl, we subtract the crystal contribution es- 
timated in Section 5.1.2 from the rate given in Table |4] For 
T3, we have no direct estimate of the crystal contribution so we 
assume conservatively that the measured rate is entirely due to 
copper surface impurities. The corresponding 90% C.L. upper 
limit for the background rate in the ROI of CUORE after the 
anti-coincidence cut is 0.02 counts/keV/kg/y for the Tl treat- 
ment and 0.03 counts/keV/kg/y for T3. If the crystal contri- 
bution in T3 could be considered to be the same of Tl, then 
the 90% C.L. upper limit for the T3 treatment would also be 
0.02 counts/keV/kg/y. These values represent for CUORE the 
major component of the background in the y8y80v-ROI. All the 
other background sources are in fact expected to give contribu- 
tions lower by about one order of magnitude (see Ref. 0). 

Conclusions 

The Three Towers Test was undertaken to validate and com- 
pare three different methods considered to control background 
due to surface radiation from copper parts in the CUORE de- 
tector. Looking for the best compromise between cost, repro- 
ducibility and background control, a mixed approach will be 
pursued for CUORE: small copper parts will be subjected to 
TECM while the large copper shields facing the crystals will be 
covered with polyethylene. These two techniques gave com- 
patible results for the event rate in the degraded a window, 
2.7-3.9 MeV, without any coincidence cuts. The data from 
the TTT apparatus demonstrate that surface contamination lev- 
els lower than 7xl0~ 8 Bq/cm 2 for 232 Th and 238 U and below 
9xlO~ 7 Bq/cm 2 for 210 Pb can be achieved. This is the best purity 
level we ever achieved for copper surfaces. An extrapolation of 



the rate measured in the degraded a window of Tl and T3 to 
the ROI of CUORE was performed. Depending on a reason- 
able range of unknown parameters in the extrapolation model 
the 90% C.L. upper limit for the ROI-background in the anti- 
coincidence spectrum of CUORE is between ~0.02 and 0.03 
counts/keV/kg/y. This is close to the design goal of the experi- 
ment and should allow CUORE to significantly improve limits 
on — or perhaps discover — Majorana neutrino masses. 
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